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Abstract: In a variant of the Randall- Sundrum braneworld model which has a charged 
black hole in an anti-de Sitter bulk, the 4D speed of gravity depends upon one's location 
in the bulk, and in general differs from the speed of light on a given 3-brane. We apply 
phenomenological constraints on the difference in the speeds of light and gravity to models 
which use the warping to solve the hierarchy problem. In particular, consideration of the 
gravi-Cerenkov radiation of Kaluza-Klein gravitons by ultra-high energy cosmic rays leads 
to a stringent limit on the fractional difference between the graviton and photon speeds in 
RS-like models: it must be less than a part in 10 23 . 
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1. Introduction 

The violation of Lorentz symmetry at high energies would be a revolutionary discovery, 
giving us important clues about the nature of physics beyond the standard model Q-Q- 
It has been argued that such violations are predicted by some theories of quantum gravity 
Q. Although string theory does not necessarily predict Lorentz violation, it can do so via 
background fields which lead to noncommutative geometry Q, or which explicitly violate 
Lorentz invariance in the gravitational sector Q. 

Nevertheless, it is not easy to invent fundamental theories which break Lorentz sym- 
metry in a plausible way. Therefore it is useful to consider new possibilities for breaking 
Lorentz symmetry. A novel approach was presented by Csaki et al. [^, in the context 
of the Randall-Sundrum (RS) solution to the hierarchy problem ||, involving two branes 
separated by an extra dimension with the geometry of 5D anti-de Sitter space. Although 
the original RS model had Lorentz symmetry, ref. [0] (hereafter called CEG) considered a 
variant in which the 5D bulk contains a black hole of charge Q and mass //, giving a metric 
(the AdS Reissner-Nordstrom black hole solution) with the line element 

r 2 1 

ds 2 = -Mr) dt 2 + -rdx 2 + TT --dr 2 ; (1.1) 
V h(r) 

Kr) = ^-^ 2 +^. (1-2) 

The extra dimension has coordinate r, and the length scale / characterizes the curvature of 
the AdS geometry, which is determined by the negative bulk cosmological constant through 
I 2 = — 6/A. (We work in units where the 5D gravitational constant = 1.) 
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The presence of a bulk black hole (BH) can be physically motivated: it was argued 
in [|i]] that emission of gravitational radiation from a brane in the early universe can 
lead to the formation of a BH in the bulk. Possible cosmological consequences of gravity 
traveling faster than light in this scenario were considered in references [11]. However in 



the cosmological solutions, the brane moves away from the BH as the universe expands, 
and as a result the presence of the BH quickly becomes irrelevant for the brane observer: its 
effects are redshifted away. The presence of charge on the BH is motivated by the desire to 
have a solution in which the brane-BH distance remains constant, and therefore its effects 
can also be important at late times. 



Along any 4D slice of constant r, the metric ( |l.l[) is Lorentz invariant, with speed of 
light given by dx/dt = (7/r)/i 1//2 . Thus an observer on a D-brane which is restricted to a 
particular value of r does not see direct evidence of Lorentz violation in the matter sector. 
However gravity propagates in all dimensions, not just along the brane. A 4D observer 
would find that gravity propagates with a different speed than matter. 

It might be thought that the weakness of gravitational interactions would render this 
mild Lorentz violation phenomenologically harmless; after all, the speed of gravity has not 
been measured. However, there are ways in which it would be manifested. One possibility 
is that self-energy diagrams with graviton loops would induce speed differences between 
observable particles, like photons and electrons (where by "speed" we mean the maximum 
speed, as the momentum becomes infinite). These kinds of constraints were considered in 

El- 

A more powerful constraint exists in the case where gravity travels slower than visible 



particles, especially protons 13]. The emission of gravitational Cerenkov radiation would 
very efficiently damp ultra high-energy cosmic rays (UHECR's) to levels below those which 
are observed. One of the main observations of the present paper is that this is precisely 
the situation in asymmetrically warped models in which we are assumed to be living on a 
negative tension brane, where the weak scale hierarchy can be solved by the warping. In 
CEG it was assumed that we live on a positive tension brane; in that the speed of gravity 
is faster than that of light, so that the present strong constraint is not relevant. 

In CEG, the asymmetric effects of warping were treated perturbatively. In the present 
work we also extend their analysis by considering these effects in the nonperturbative 
regime, and for the Kaluza-Klein (KK) excitations of the graviton. These could be relevant 
sources of Lorentz violation in colliders, since in RS-like models the KK gravitons couple 



strongly to the standard model particles [14]. 



We begin by reviewing the CEG model and their perturbative derivation of the speed of 
gravity, and we enumerate the most stringent experimental constraints on Lorentz violation 
in these models. In section 3 we point out that the perturbative treatment can be misleading 
if the observer is located on the positive tension brane, whereas it gives reliable predictions 
for observers on the negative tension brane. In section 4 we show that the KK excitations 
of the graviton have similar Lorentz-violating properties to the zero mode, considerably 
strengthening the Cerenkov bound introduced in section 2. 
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2. The Models; perturbative speed of gravity 



2.1 Model with two branes 

Similarly to the RS model, CEG cuts the bulk at two positions, call them r + and r_, by 
the insertion of branes and the use of Z2 orbifold symmetry to define the discontinuity in 
the derivatives of the metric functions at the branes. The parts of the solution r > r + and 
r < r_ are discarded, and replaced by mirror copies of the kept part of the solution to create 
a solution with orbifold fixed points at r±. Unlike RS, it is not possible to use conventional 
branes with equation of state p = —p at both positions. There are four junction conditions 
which relate the energy densities p± and equations of state p± = w±p± to the black hole 
parameters: 

pi 2 = 3 (l + |±p^ 2 ) r 4 = 3 (l + ^ P 2 J 2 )ri; (2.1) 

12,2 _ o ft , l + 3w + 2 ,2\ R n /1 j 1 + 3U>_ 2 2 \ 6 



Q'Z' = 2 \1 + ——L p ^y* + = 2 I H -^—pt.l \r_. (2.2) 

It is easy to see that taking w+ = W- = — 1 implies that r + = r_. Let us consider the 
physically interesting case of a large hierarchy, where 

e = r_/r+ - 1(T 16 . (2.3) 

Physical masses on the negative tension brane at r_ (the "TeV brane") are suppressed 
relative to their bare values by the factor e, giving a resolution of the weak scale hierarchy 
problem ||. By solving the jump conditions ( |2.1| )-( |2T2] ) for e, we see that there are two 
ways of achieving a small value of e: either by (1) tuning the numerators (1 + ^-p\l 2 ), 
(l + ^g^/?2j 2 ) to be small, or (2) tuning the denominators (1 + ^p 2 J 2 ), (l + ^^p^ 2 ) 
to be large. 

1. In the first case, lp-/Q and W-(lp-/Q) 2 can naturally be of order unity, while u>+ and 
p + must be tuned to the values 

A , w - 2 ,2 



w+ ^_l_ 2e 4 (l + — ptl 1 ) <-l; Zp+=6 + 0(e 4 ), (2.4) 



implying the relation 



Q 2 Z(l + ^p 2 J 2 ) 



p = Cl \; c x = > J> ' > = 0(1). (2.5) 
2(1 + ^^) 

The inequality in (^4|) is a violation of the weak energy condition (albeit a small 
one), which states that p > —p, and would thus require some exotic kind of stress 
energy on the Planck brane. 

2. In the second case, we need 

w_ = 0(6 2 ); /p_ = 0( e - 6 ), (2.6) 

which implies 



91. Cz _ 3(l + ^ 2 

r 2 ' 2 2 (l + i+g± p ^2 



P = C2^; c 2 = "J . \ =Q(i). (2.7) 
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The relations (|2,5| ) , ( p. 7] ) are useful for computing the speed of electromagnetic radiation 
(or any relativistic particles) on the branes. Working to first order in the Lorentz violating 
parameters fi and Q 2 , the deviation in the speed relative to unity is 

Hi* , Q 2 l 2 PQ 2 / 1 *\ 

where the choice ± refers to which brane one is on, and the label i refers to the choice of 
tunings immediately above; hence r\ = r_ and ri = r + . 

2.2 Model with one brane and a horizon 

The preceding discussion assumed the existence of two branes bounding the extra dimen- 
sional space. It is also possible to eliminate the TeV brane altogether by not cutting out 
the small r region containing the black hole. In that case one should choose \i and Q 2 
so that h(r) vanishes at least once between r = and r = r + . This insures there is a 
horizon shielding the black hole, in accordance with the cosmic censorship hypothesis. In 
this scenario we would necessarily be living on the Planck brane at r = r + . The parts 
of the jump conditions (^1])-( |2^ ) involving p + and w + would still be valid. CEG showed 
that it is not possible to embed the brane and still have a horizon shielding the black hole 
unless the weak energy condition is again violated, w + < — 1. This is most easily studied 
when the special relation 4/i 3 / 2 = 27Q 4 holds; this is the limiting case where the two hori- 
zons of the AdS Reissner-Nordstrom black hole degenerate into a single horizon, located at 
r h/l = (2Q/Z 2 ) 1 / 3 /v / 2- One can solve the jump conditions to find that the brane's equation 
of state depends on its position, r+, and is given by 

w+ =-i- 2 (^y +o(^y . (2.9) 

Hence the violation of the weak energy condition becomes small as the brane is moved 
farther away from the black hole. Ref. [jl^] showed that this violation could be moved into 
the bulk stress-energy tensor, but never eliminated, so long as the horizon exists. The 
deviation in the speed of light on the brane in this scenario is given by 



fc » = -H^J H^J • (2 - 10) 

We mention this case for completeness, but we will not further consider how gravity prop- 
agates in this case. 

2.3 The speed of gravity; experimental constraint 

Now we turn to the speed of gravity. Its determination is simplified by CEG's observation 
that the dynamics of gravity are the same as those for a massless bulk scalar field, with 
action and equation of motion 

S = 1 f d A xdry/^g AB dA4>d B (p] (2.11) 



2 

SS , — w 



d A (V=99 A dB(j>) = 0. (2.12) 
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u = \q\+Su ; 5a;o = ^(-2^ + 5 r (l + e- 2 )). (2.13) 



The problem can be further simplified by splitting the metric into a pure AdS part and the 
perturbation involving the black hole parameters. In other words, one works perturbatively 
to first order in [i and Q 2 . The solution for the nth KK mode is separable, (p n (x^,r) = 
e %<yUJnt ~^' x ^ (f) n {r) , where 4> n ( r ) satisfies Neumann boundary conditions at the branes, since 
the jump conditions have already been satisfied by the background solution. The 

dispersion relation of the graviton zero mode solution in the tower of KK excitations can 
be determined analytically (in the model with two branes) for a solution with 3-momentum 

Mi! (_ojL + 9. 2 

4e 2 ^ r 4 t- r e 
The deviation in the speed of gravity from 1 is given by 

= |f- (2-14) 

The relevant Lorentz- violating observable is the difference 5c em — 5c grscv . In general, the 
sign of this difference depends on the value of the 0(1) constant c\ or C2- However for TeV 
brane observers, if the hierarchy is very large (e <C 1), and if we exlude the case (|2.4|)-( |2lj| ) 
which violates the weak energy condition, the difference is completely dominated by the 
term / 2 Q 2 /(2r 6 ) in (pT8|) , and leads to 



l 2 Q 2 

2r 6 • 



Because the difference is positive, the stringent constraints due to gravitational Cerenkov 
radiation of UHECR protons is applicable [0| : 



Sc em - 5c grav < 2 x 1(T 15 (1(T 19 ). (2.16) 

where the less stringent limit assumes the UHECR's are of galactic origin, and the more 
stringent one applies if they originally come from neutrinos originating from cosmological 
distances. Therefore we have the bound 

6 -^4 - ^ < 4 x 10- 15 (10- 19 ). (2.17) 
n(r_) r_ 

Notice that the quantity on the left is the fractional perturbation to the metric function 
h{r) due to the dominant Lorentz violating term, evaluated on the TeV brane: h(r) = 
(r/l) 2 (l + 5h/h) It is therefore natural that this is the combination which is experimentally 
bounded. (l2~17|) is one of the new constraints derived in this paper. 

If we are willing to entertain the possibility of case 1, with exotic matter on the Planck 
brane which violates the weak energy condition, then the difference of speeds as measured 
on the TeV brane becomes dependent on the details of the TeV brane stress energy, 



/2Q2 

(5c em - ^c grav = ^-g-(l - ci). (2.18) 
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If ci < 1, the Cerenkov bound again applies, in the obvious way. If c\ > 1, gravitational 
Cerenkov radiation cannot be produced. Then the strongest bounds come from tests for 
deviations from general relativity in the parametrized post-Newtonian (PPN) formalism. 
A difference between the speed of light and that of gravity implies the existence of a 
preferred reference frame, the one in which we have been working, where the speeds of 
gravity and light are independent of direction. In a boosted frame, one of these would 
become anisotropic, depending upon whether the boost respects Lorentz invariance of the 
gravitational or the electromagnetic propagation. (One but not both can be preserved, due 
to the difference in speeds.) In such a situation, the parameter oli of the PPN formalism 
is nonzero. Essentially, we have two metrics, one for photons and one for gravitons, so the 



effective theory is Rosen's bimetric one [16]. Such theories predict a torque which would 
cause precession of the sun's spin axis. The latter is closely aligned with the solar system's 
planetary angular momentum vector. If we assume that the alignment is not coincidental, 
then precession due to PPN effects is constrained, and leads to the bound [17] 

|aa| =2 |fc cm -fc grav | < 12 x 1Q _^ (2 ig) 

which carries over to the quantity —$-11 — ci|. 



3. Nonperturbative analysis 

The above discussion assumed that the dispersion relation of the graviton gets modified in 
the simple way which was predicted by treating the Lorentz- violating part of the metric, 5h, 
to first order in perturbation theory. However, at large graviton momenta, the perturbative 
prediction can be modified. To see this, one must examine the equation of motion ( 2.12] ) . 



Let us rewrite it to first order in Sh, with <ft^ and respectively denoting the zeroth 
and first order solutions, in powers of Sh, for the Kaluza-Klein zero mode. We will also 
change coordinates to the form r/r+ = e~ ky with k = The wave equation for the zero 
mode with 3-momentum q is 

& v $\y) - Ak sgn(y) d y $\y) = e 2k ^ (^qu™ + q 2 ^ ( 3 .1) 

= e 2 %l (-2,4 1) + fJ(-^ 4%l + ^%l))4° ) . 



Regardless of how small Sh/h is, at sufficiently large momenta the source term q 2 (5h/h) < 
becomes so large that it can no longer be reliably treated as a perturbation. 

To explore this, we have numerically solved the full, unperturbed graviton wave equa- 
tion to find the dispersion relation for the zero mode, uo(q), and compared this result to 
the perturbative prediction ( |2.13| ). To make the problem tractable, we have chosen definite 
values for the brane stress energy parameters. For the tuning of parameters, we adopt 
case 2 above, eqs. Q2.6[) - (^?7|) , so that C2 depends on the positive tension brane parameters. 



1 We thank Guy Moore for discussions on this point. 
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These we take to be w+ = — 1, while leaving p + free to vary. The full wave equation can 
be written as 



(i s p 6k\y\ \ 
t,\ ^)+ e2% ' 
h(y> J 



U! 



h 2 (y) h(y) 



<t>(y) = o, 



(3.2) 



where 



h(y) = l + 5e Ak \ y \ (-2, + 2e 2k \A ; 5 = 1 



36 



(3.3) 



In this form, it is clear that the Lorentz-violating effects, coming from h— 1, are maximized 
at the TeV brane, where y = ?/_. For the purposes of comparing exact results to the 
perturbative ones, we quantify the perturbation by the parameter 



A = 25e 6ky ~ . 



(3.4) 



We expect the perturbative treatment to be valid when A <C 1. 

With these definitions we can now compare the exact numerical results with the per- 
turbative ones. The deviation of the dispersion relation of the zero mode, relative to its 
Lorentz-conserving value, is plotted as a function of momentum for several values of A, 
and as a function of A for several momenta, in figure |l], where Su> and q are given in units 
of k = We have taken a very small hierarchy, ky_ = 2 for this illustration, and will 
comment on the extrapolation to larger values below. The interesting feature is that the 
true 5uj(q) flattens to a constant value at large q, contrary to the perturbative expectation. 



q 

Legend 



A=.10 
A=. 30 
A=. 70 
A=l .00 



A 

Legend 



q=.05 
q=. 65 
q=1.35 
q=2 . 00 



Figure 1: Correction to uj^ /k as a function of q/k (left) and A (right). The upper (straight) lines 



show the perturbative prediction (2.13), while the lower (curved) lines show the numerical result. 



Each color represents a different perturbation or momentum, and y_ = 2/k. 

If we compare the speed of gravity to that of hypothetical photons which are trapped on 
the positive tension brane, a consequence of this behavior of the graviton dispersion relation 
is that the speed difference (with speed defined using the group velocity v = duj/dq) does 
not remain constant at large momenta, but rather vanishes as q — > oo, as shown in figure 
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IH At large momenta, gravitons tend toward the same speed as radiation on the Planck 
brane. This can be understood as follows. The graviton zero mode is localized on the 

2 2 

Planck brane. At large q 2 , the effect of the terms 4r) m the equation of motion 

h 2 (y) Kv) 

Q3.2| ) is to localize it even more, thus driving the graviton to resemble more closely radiation 
which is trapped on the Planck brane. This trend becomes evident for momenta q > 1/7, 
where the speed difference is significantly reduced relative to its maximum value at q = 0. 

As the hierarchy between the Planck brane and the negative tension brane is increased 
by letting y_ go to larger values, the Lorentz violating effects seen by a Planck brane 
observer are suppressed, if the parameter A is held fixed. The magnitude of c grav — c em 
scales like e~ 2ky ~ . This can be understood analytically, using the perturbative results of 
the previous section. One can show that, for the present choice of parameters, 

c gr av -c em = ^(cosh(2%_)-l)e- 4 ^-A ~ +0{e' 2k y-). (3.5) 

At the same time, the g-axis in figure || gets rescaled by a factor of e~ ky ~ , such that 
the range of momenta with an appreciable deviation in c grav — c em gets shifted to smaller 
physical values. In the hierarchy-solving case e~ ky ~ ~ 10 16 , this corresponds to the TeV 
scale. 




Figure 2: Relative speed difference between the gravity zero mode and a photon on the positive 
tension brane as a function of q (left) and A (right) when the physical brane has positive ten- 



sion. The upper (straight) lines show the perturbative prediction derived from (2.14), while the 
lower (curved) lines show the numerical result. Each color represents a different perturbation or 
momentum, and y_ = 2/k. In the right panel, all the perturbative curves coincide with each other. 

On the other hand, an observer on the TeV brane will not be concerned with this small 
momentum-dependence in the graviton speed, because the difference between c grav and c em 
remains relatively large even at high momenta, as shown in figure ||. Moreover, in contrast 
to the case of the Planck brane observer, the difference c grav — c em remains constant as the 
hierarchy is increased, if A is held fixed. Again, this can be understood by computing the 
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speed difference perturbatively, for the present choice of parameters, 

Cg r av-Ce m = - (cosh(2%_ ) - l)e~ 2ky ~ A 0(1) x A. (3.6) 



-0 . 12 -j ,. . — 

-0 .16 
-0 . 13 

-0 .2 : 
-0.22- 
-0 .21- 
-0 .26: 

-0.2E- 

0.5 1 1.5 2 

Legend ^ 
A=.30 

A=. 50 

A=.80 
A=1.00 



Figure 3: Same as figure 2, but here the photon is on the negative tension (TeV) brane. All curves 
are essentially coincident in the right panel. 

The conclusion of this section's analysis is that we can trust the results of the per- 
turbative treatment if we assume that observers are living on the negative tension brane, 
as one would expect if the hierarchy problem is being addressed. Only for Planck brane 
observers would it be important to distinguish the perturbative from the exact results. 
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4. Lorentz-violating KK modes 

It would be interesting if Lorentz violating kinematics could be observed in the laboratory, 
in particle collider experiments. In section 2 we noted that the speed difference between 
gravity and light on the TeV brane should be less than a few parts in 10 7 . Such a small 
effect would probably be impossible to see at the LHC, even though KK gravitons inter- 
act strongly (with TeV-suppressed instead of Planck-suppressed couplings) with standard 
model particles and would be copiously produced, if sufficiently light. 

Even if their Lorentz-violating properties cannot be directly detected, an interesting 
indirect effect is the Cerenkov emission of KK gravitons from UHECR's, which would 
strengthen the bound mentioned in section 2. That bound conservatively counted only the 
damping due to emission of massless gravitons. But KK modes can also be emitted, so 
long as their mass obeys the inequality 

2 

< c em - c grav (p), (4.1) 

where p is the momentum of the UHECR. This is the condition for Cerenkov emission to 
be kinematically allowed. In the situation where the warp factor e = 10 -16 , the KK mass 
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gap is of order TeV. The highest energy cosmic ray which has been detected had energy 
3 x 10 11 GeV [18]; thus for a speed difference that saturates the solar system bound ( 2.1SQ , 
the number of relevant KK modes is of order 10 5 . It is therefore worth exploring whether 
the KK modes have similar Lorentz-violating kinematics as the graviton zero mode. 

We have addressed this problem both numerically and analytically. The analytic ap- 
proach is to solve the wave equation once again treating 5h to first order in perturba- 
tion theory, but now expanding around the zeroth order Lorentz-conserving KK wave 
function, (j) n {y) = 4>n \V) + 4>n\y), in order to find the corresponding energy eigenvalue, 
uj n (q) = LUn\q)+LL>n\q). The solution is described in detail in the Appendix. We have also 
numerically solved the full wave equation using the shooting method, to find the dispersion 
relation uj n (q). The deviation of u) n (q) from the standard result u = q is shown for the 
first two KK modes in figure ^ The comparison between the perturbative and numerical 
results is shown there. 



o.sj 



I D ) 
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0.5 1 1.5 
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Legend 

A=.2 4 
A=. 4 8 
A=. 7 6 
— A=1.0 



Legend 
q=0 
q=. 95 
q=2 .00 
q=3.Q0 



Figure 4: Same as figure 1, but for the first two massive KK modes. The perturbative prediction 
comes from eq. ( |A.8| ). 

As we did for the KK zero mode, we can compute the difference in the speed of gravity 
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relative to that of particles trapped on the positive tension brane. In order to highlight the 
differences which are due to Lorentz violation, we take the conventional trapped particles 
to have the same rest mass as that of the KK graviton to which it is being compared. 
The results, shown in figure |5] are qualitatively similar to those for the zero mode: gravity 
is faster than particles on the Planck brane, but the speed difference becomes smaller at 
higher momenta. 




Figure 5: Relative speed difference between the first two massive gravity KK modes and hypo- 
thetical particles of the same masses, living on the positive tension brane, as a function of q (left) 
and A (right). Upper curves are perturbative, lower are numerical results. 



More interesting is the speed difference relative to particles on the TeV brane. Figure 
|6] shows that, again like the zero mode, KK gravitons are slower than same-mass particles 
on the TeV brane. Therefore our expectation that they are produced in gravi-Cerenkov 
radiation by particles with momenta satisfying (4.1) is justified, and we should revise the 
bound. Adapting the result of [13] to the tower of TeV-mass-gap gravitons whose couplings 
are only TeV-suppressed, the rate of energy loss of a UHECR with momentum p and mass 
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Figure 6: Relative speed difference between the gravity KK modes and hypothetical particles of 
the same masses, living on the negative tension brane, as a function of q (left) and A (right). Upper 
curves are perturbative result, lower are numerical. 



-7^773/ dgq I dM sm±e, (4.2) 



m p is given by 

ilX 

tft TeV 3 7o Jo 

where g is the momentum of the emitted KK graviton, M is its mass, and 6 is the angle 
of emission of the Cerenkov radiation. The kinematic limits and 9 (in the limit 9 <C 1) are 
given by 

9 2 = (l-q/p)q- 2 (Ml ax -M 2 ) (4.3) 
ML* = Q 2 25c- 3- L (4.4) 



p 2 (l - <?/p) 
2p 2 5c 



V U-T^k ■ (4-5) 
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Here 5c = c em — c grav , which must be positive in order for Cerenkov emission to occur. 
Evaluating the integrals gives 

^ ~ 0.1 (5cf' 2 f{5c) (4.6) 
dt \ j j\ / Te y 3 v > 

where the function / is approximately 1 for 5c 3> ^ and drops very rapidly to zero (like 

2 

(1 — m 2 /2p 2 '5c) 11 / 2 ) as 5c — * The bound on 5c comes from demanding that dE/dt be 
less than p/L for a cosmic ray which is propagating over a distance L. For the highest 
energy cosmic ray observed, which is identified with a proton, the bound is satisfied by 
taking 5c to be below the kinematic limit where / = 0. This implies that 

2 

5c < ^ 9* 1(T 23 . (4.7) 
2jr 



5. Conclusions 



Asymmetric warping can provide a plausible means of introducing Lorentz violation into a 
theory with extra dimensions, which is essentially a form of spontaneous breaking due to the 
gravitational background. Since at tree level the violation is confined to the gravitational 
sector, the effects can be sufficiently weak to be at the borderline of detection. In this paper 
we have explored some of the consequences of a theory where the asymmetric warping comes 
from the presence of a charged black hole in the extra dimension. Such a scenario can be 
compatible with the RS solution to the weak scale hierarchy problem. Unlike RS, in the 
black hole case it is necessary to allow for nonstandard equations of state for the tensions 
of at least one of the branes, though it is still possible to respect the weak energy condition. 
One issue we have not explored is the role of a stabilization mechanism such as that of 
Goldberger and Wise |l9| for the size of the extra dimension. One might expect that the 
bulk is already stabilized once the brane energies and equations of state are fixed, since 
algebraically the ratio of brane positions r_/r + is determined. However it is possible the 
radion mass 2 is negative, even though it is not zero. Even if it has the right sign, the 
distortions of the bulk geometry by a bulk scalar could have an effect on the propagation 
of gravity. It would also be interesting to know whether pure tension branes could be 
admitted with the addition of a bulk scalar, as it would be desirable to eliminate the need 
for an unusual equation of state. Finally, we did not consider how Lorentz violation is 
manifested in the single brane case, mentioned in section 2.2, although we expect it to be 
similar to the case of a Planck-brane observer in the two brane model. These subjects 
could merit further study. 
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A. Perturbative KK mode solution 



The zeroth order in Lorentz- violation KK graviton wave function (which like the zero mode 
obeys Neumann boundary conditions) is pi] 
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(A.l) 



where N n is a normalization constant. The mi are determined by the boundary conditions; 
in the limit of a large hierarchy, they satisfy J\(m n e ky - /k) = 0. 

The wave function correction 4>n (y) is difficult to compute, but as is familiar from 
perturbation theory in quantum mechanics, we don't really need it if we are only inter- 
ested in the correction to the eigenvalue, ujn(q)- In the quantum mechanical analogy, 
we require only the matrix element of the perturbation to the Hamiltonian. Defining 
V'n(y) = e _2fc ! y !</> n (y), the equation of motion becomes 



HMv) 



where, to first order in the perturbation, 



V'n(y), 



H = dy(h(y)dy) - 4k 2 h(y) + 2k sgn(y) h'(y) - q 2 e 2k ^ + 4kh(y)d(y) 
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(A.2) 

(A.3) 

(A.4) 
(A.5) 



with =h-l. 

We take the inner product with ip^ and use the Hermicity of to cancel the terms 
involving ipn\ This yields 

{1) _ -rf^W|4 0) ) + (^ ) ) a (4 0, |fc (1) e a *"»W ) ) 



2u$>(^\e*M\$>) 



(A.6) 



The inner product in the denominator is unity, with the appropriately normalized wave 
function. Putting in the explicit form of H^> , the expression becomes 



(i) 



C (0) (dyCh {1) dy) - 4k 2 h^ + 2kh'W) 4 0) dy 



+4°) / ~ h^e 2k y\^\ 2 dy 



.— A + J^B 
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(A.7) 
(A.8) 



Since iOn^ = \/ m n + q 2 J and m n /k oc e ky ~ , all the terms are roughly comparable for 
q <C m n . But, for large momentum q, only the second integral is important. In that case, 

ry- 



(1) 



\^\ 2 h {1) e 2ky dy 



(A.9) 
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k 2 m 



jl~m 2 [-^ + ^y y (A.10) 

^;<^)(-r y+ f e i^ (A - n) 



where the last approximation holds in the limit of a large hierarchy. The resulting integral 
can only be done numerically. However, our analysis is still useful since it tells us that the 
first correction term to u n is proportional to q for large momentum, just as for the zero 
mode. Moreover we can show how the dispersion relation is expected to change. Squaring 
u n gives 

u? n = (4°) + u>W + • • -f = (ml + q 2 )(l + 2B n - 2A n ) (A.12) 

to first order in the perturbation. Comparing with the usual dispersion relation, one sees 
that the limiting speed of the nth mode has been modified to 

4av n = l + 2#n, (A.13) 

and that the mass of the mode changes due to the presence of the black hole perturbation: 

Ml = m 2 n (l + 2B n ) - 2A n . (A.14) 
In terms of these, the group velocity of the nth mode with momentum q is 

%av n = ^ ° l ^ q ■ (A.15) 

As a consistency check, it can be verified that this result reduces to the previous one 
for the zero mode case where it is possible to do the integrals. 
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